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ABSTRACT 

Inspired by the results of the Ising model within an imaginary external 
magnetic field, we introduce a transformation in quantum systems with a 
^-vacuum term that amounts to a rescaling of z = cos(|). Making use of this 
transformation we are able to determine the order parameter as a function 
of 9. The approach is successfully tested in models with both broken and 
unbroken CP symmetry at = it. 
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Quantum systems with complex euclidean actions have received an im- 
portant amount of attention in recent time. Indeed there are two very rele- 
vant subjects in quantum field theory the analysis of which requires the use 
of complex actions: the study of topological structures and the behavior of 
matter at high temperature and density. Notwithstanding the big interest 
in these subjects, it is currently not possible to write an efficient code to 
perform numerical simulations of systems with complex actions. All efficient 
algorithms to simulate statistical systems are based in the use of a true prob- 
ability distribution function as Boltzmann weight. This is the reason why all 
attempts to analyze systems with complex actions are based (i) on the use of 
specific properties of the particular model under study [T], (ii) uncontrolled 
analytical continuations, or (Hi) in a very delicate process of summing up all 
the contributions which appear in the definition of the partition function [2] , 

0, 0, 0, 0. 

In a recent paper [JJ we proposed a new approach to analyze systems with 
a ^-vacuum term, tested it in two analytically solvable models, and applied 
to the study of CP 3 with a topological term in the strong coupling region. 
This approach is based on the use of a high accuracy method to measure the 
probability distribution function of the topological charge density at 9 = 
from numerical simulations at imaginary 9. It is combined with a multipreci- 
sion algorithm which allows to sum up all the terms which appear in the 
partition function, that differ by many orders of magnitude, running from 1 
to e~ v , where V is the lattice volume. 

Even if the results reported in [7j are quite promising, the method is not 
free from systematics and it would be therefore important to have alternative 
approaches available in order to compare physical predictions. This is indeed 
the purpose of present paper. 

We should say from the beginning that contrary to the approach reported 
in [JJ where no assumption on the phase structure of the system was made, 
here we will make the assumption that the system has no phase transition at 
real 9 except at most at 9 = n. This means that our proposal should be used 
in practice as complementary to other approaches and one should look for 
consistency between their results. 

The method is based in extrapolating a suitably defined function to the 
origin. The function turns out to be very smooth in all the cases we have 
considered up to now and this makes us confident on the whole procedure. At 
the very end it is equivalent to an analytical continuation from imaginary to 
real 9. However it has not to be confused with analytical continuation of the 
free energy density and/or the order parameter which are extremely sensitive 
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to statistical errors and to the choice of fitting functions and therefore with 
small or no practical utility. 

To fix the ideas and to keep this paper selfcontained, let us write the 
general form of the partition function of a system with a topological ^-term 
as a sum over all topological sectors of the partition functions for each sector 
weighted with the proper topological weight, 

Zv{0) = J2Mn)j 6n , (1) 

n 

where py(rz) is, up to a normalization constant, the probability of the topo- 
logical sector n at 6 = 0, and the sum runs over all integers n. Writing 
equation (JTJ) one needs to assume that CP symmetry is realized at 6 = 
since otherwise Z v {6) would be ill-defined jH]- 

In almost all practical cases the sum in has a number of terms of order 
V since the maximum value of the topological charge at finite volume is of 
this order. The partition function (0) can also be written as a sum over the 
density of topological charge x n = n/V, and setting pvij 1 ) — exp[— V fy(x n )] 
one gets: 

Z v {6) = e- yfv{Xn) e WVx " . (2) 

Equation (j2J) defines a 2tt periodic function of 9. For models CP-invariant 
at 6 = fy(x n ) is an even function; this property allow us to write: 

cos 2 -j (3) 

yn 

where y n , which is non negative by construction, takes the same non negative 
values as x n and Gv{y n ) is a functional of fv{x n ). The mean value |/(cos|) 
of y n computed from (jHJ) is related to the mean value x{6) of i times the 
density of topological charge as follows: 

cos - j tan - (4) 

which shows explicitly the 2tt periodicity. At imaginary 6,6 = —ih, the 
previous equation becomes 

x(—ih) = y ^cosh — j tanh ^ (5) 
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Furthermore, at real h, the system described by (JTJ) can easily be simulated 
by standard numerical algorithms since in this case the Boltzmann weight is 
well defined. From numerical simulations at real h one can get y(z) (jSJ) for 
values of the argument z = cosh | between 1 and oo. In order to get the 
topological charge density at real 6 one should extend this function to the 
(0, 1) interval in z. But as stated before, an extension of y(z) is sensitive 
to statistical errors and fitting functions so that it makes this procedure not 
reliable. 

The function y(z) depends also on the couplings of the model under con- 
sideration. For the one- dimensional Ising model within an external real or 
imaginary magnetic field h, the function y(z) can be analytically computed 
and the final result is 

V& = j * 1/2 (6) 

(1 + (e-^ - I) -1 z 2 ) 1 

which shows explicitly how this model, at imaginary values of the magnetic 
field (0 < z < 1) and antiferromagnetic couplings (F = < 0) breaks 
spontaneously the Z{2) symmetry at h = in (see (JH)). But this model shows 
also a very interesting property: a change in the coupling F, or what is the 
same, a change in the physical temperature T, is equivalent to a rescaling 
of the variable z since y depends on z through the combination e^z, with 
A = — log(e _4F — 1). In other words, if we define the transformation: 

y x (z) = y(eh^j (7) 

which would allow us, for negative values of A, to extend the order parameter 
from imaginary 9 to real 6, the 1 — d Ising model shows the exceptional 
feature that the transformation (|7J) is equivalent to a change of the coupling. 
This property, if true for other models, would be extremely interesting since 
it would reduce the problem of complex actions to simulations with real 
actions. Unfortunately this seems not to be the case since it is easy to check 
that already compact U(l) in two dimensions with a #-term does not exhibit 
such a property. However this is not the end of the story and we will try in 
what follows to make use of transformation (J7|) in a much less ambitious way. 

Many of the physically relevant models, as CP^ -1 or QCD, are asymptot- 
ically free. Asymptotic freedom implies that the continuum limit is reached 
in the weak coupling region, but in this region the order parameter (and 
y{z)) will be very small for values of z near 1 since the density of topological 
structures in lattice units is very suppressed at weak coupling. This means 
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that, from simulations at real h, we can plot y\/y against y up to values 
of y very close to the origin, provided that y(z) does not vanish for z > 0. 
We cannot exclude this possibility completely, but it does not happen in 
any of the exactly solvable models we know. Indeed, it is easy to show that 
y(z) 7^ for z > if all the coefficients Gv(y n ) entering the right hand side of 
equation (j2J) are positive, as it happens in the models studied in the present 
work. In such the fluctuation-dissipation theorem ensures that y(z) 

is a non decreasing function of z. Furthermore, if the order parameter does 
not diverge at 9 = tt, y(z) should vanish at least linearly with z at z = 
(that corresponds to 9 = tt). 

Moreover, if y\/y as a function of y shows a smooth behavior near y = 0, 
we can expect a simple extrapolation to y = being reliable. This point is 
crucial in what follows and obviously can not be taken for granted. Apart 
from being reasonable it has been verified in all the models whose analytical 
solution is known to us. In particular it is true for the one-dimensional 
Ising model, compact U(l) in two dimensions, CP 3 , and a toy model that 
exhibits no phase transition at 9 = tt, thus covering all the possible different 
qualitative possibilities. 

In the Ising case we have done simulations of the model within an external 
real magnetic field h at F = —2.0. This coupling was chosen in order to be 
in the "weak coupling region" where y(z) is small near z = 1. The model 
breaks spontaneously the Z(2) symmetry at 9 = —ih = tt and therefore y(z) 
vanishes linearly with z (see (jlj, (JJJJ)). 

If we define an effective exponent 7a by the following equation: 



this exponent, in the limit y(z) — > (z — > 0), will give the dominant power 
of y(z) as a function of z near z = 0. A value for this exponent that goes 
to 1 when y — > will imply spontaneous symmetry breaking, and values 
between 1 and 2 will signal a second order phase transition at 9 = tt with 
a divergent susceptibility. The reason is that the order parameter, x = 
tan(0/2)y[cos(0/2)], will behave as (tt - 6)^~ l for 9 -> tt. 

Fig. 1 reports y\/y against y for the Ising model at A = 0.5. The solid 
line is a quadratic fit of the points. This figure contains also the results for 
the effective exponent 7a against y (lower points). A simple quadratic fit 
predicts 7a(0) = 0.997(3) i.e., spontaneous symmetry breaking. We have 
checked that, as expected, we get compatible results using different values of 
A = 0(1), and the same is true for the rest of the models considered. 




(8) 
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The fit of yx/y against y together with the dependence of y(z) on z 
for z > 1 (computed from numerical simulations at real h) allows one to 
reconstruct the order parameter at real 9. To simplify the discussion let 
us assume we have computed y(zo) = y°. Next using the fit of Fig. 1 we plot 
yx against y and look in this plot for the value of y = y l which corresponds 
to yx = y° (see Fig. 2). It is simple to verify that y 1 = y(e~ x / 2 Zo). The next 
step is to iterate this procedure: we replace y° by y 1 and get y 2 = y(e~ x z ), 
and so on. Once we have computed in this way y(z) for as many values of 
z < 1 as desired, the order parameter can be computed using equation (jlj). 

Fig. 3 contains the results for the Ising order parameter. We report 
also in this figure the exact analytical results in order to put in evidence the 
agreement. 

For two-dimensional compact U(l) with topological charge we chose to 
work at /3 = 1.0 and A = 0.5. Figures similar to 1 and 3 were obtained for 
this model. The agreement between numerical and analytical results for the 
order parameter was even better than for the Ising model. Because of space 
limitations the corresponding figures are not included. 

The next model we want to analyze here is CP 3 with topological charge. 
Contrary to the Ising and two-dimensional compact U(l) cases, in which we 
can compare with the exact analytical results, CP^" 1 models, which share 
many interesting properties with SU(N), are not analytically solvable. 

We have adopted for the action the "auxiliary U(l) field" formulation: 

Zn+iiZnUn,ii ^n^n+/i^n,/j 2) (9) 

n,/i 

where iV-component complex scalar field that satisfies z n z n = 1 and 

U Utfl is a U(l) "gauge field". The topological charge operator can be con- 
structed directly from the U(l) field: 

p 

where U p is the product of the U(l) field around the plaquette and — n < 

l0g(£/p) < 7T. 

Simulations were done at f3 = 0.4. This value of the coupling is not in the 
weak-coupling scaling region, but the values of y(z) near z = 1 are smooth 
enough to allow a simple extension of yx/y until y = 0. Fig. 4 contains the 
effective exponent 7a against y at A = 0.5. A simple quadratic fit predicts 
7(0) = 1.03(5) i.e., spontaneous symmetry breaking at 9 = it. By imposing, 
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as suggested by the previous fit, 7(0) = 1.0 in the fit of y\/y against y, we get 
for the order parameter the results reported in Fig. 5. The continuous line 
in this figure represents the numerical results obtained within the approach 
reported in [7j. The very good agreement between the two approaches makes 
quite reliable the numerical results and CP structure at 9 = tt of Fig. 5. 

The models previously analyzed have a common feature: CP symmetry 
is always spontaneously broken at 9 = tt. In order to check that our approach 
works also in the symmetric case we have analyzed a simple toy model, with 
only one effective degree of freedom, where CP symmetry is enforced at 
9 = tt. This model, which resembles very much the free instanton gas model, 
is described at finite volume V by the following partition function: 



The order parameter at real 9 



shows explicitly CP symmetry at 9 = tt, as it should be since this is a model 
of one degree of freedom. We have not done for obvious reasons numerical 
simulations of this model at real h = —i9 but in order to work in realistic 
conditions, we modified the exact value of the order parameter at real h 
by adding to it a random gaussian error of 1 percent. Fig. 6 contains the 
exact analytical results for the order parameter at A = 0.005 and the points 
obtained within our approach, using always for the parameter A defined in 
© the value 0.5. The agreement is quite good. 

A final comment on the possibility of applying this approach to the inverse 
problem i.e., finite density QCD, is worthwhile. QCD at imaginary chemical 
potential and staggered fermions can be simulated numerically because the 
fermion determinant is positive definite. Its partition function belongs to the 
general class of partition functions described by (J2J), where x n now stands for 
the density of baryonic charge and exp[— Vfy(x n )] is the probability distri- 
bution function of the baryon density at vanishing chemical potential. From 
numerical simulation of QCD at imaginary chemical potential one can mea- 
sure the mean value of the density of baryonic charge and from it we can get 
y(z) for z- values between and 1. In order to get the density of baryonic 
charge at real chemical potential we should extend y(z) to ^-values larger 
than 1. Summarizing we know y\/y as a function of y around y = and we 
need to extend this quantity to larger values of y. But since a phase transi- 
tion is expected in QCD at finite chemical potential, a singular behavior of 
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U\/y as function of y should also be expected, and any analytical extension 
of this quantity should fail when trying to reproduce the correct behavior. 

This work has been partially supported by an INFN-CICyT collaboration 
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Fig. 1 
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Figure 1: y\/y and 7a as a function of y for A = 0.5 in the Ising model 
(F = -2.0). 
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Fig. 2 
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Figure 2: An example of the iterative procedure to reconstruct y(z): y\ as 
function of y is plotted using dashed line in the region where direct mea- 
surements are available and continuous line in the extrapolated region. The 
straight continuous line is y\ = y. 
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Figure 3: The order parameter in the Ising model (F = —2.0). Continuous 
line is the exact result, points have been obtained from quadratic fit of y\/y 
data reported in Fig. 1. 
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Fig. 4 
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Figure 4: CP 3 model at (3 = 0.4 in a 100x100 lattice: the effective exponent 
7a versus y. The continuous curve is a quadratic fit to the data. 
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Fig. 5 
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Figure 5: CP 3 model at (3 = 0.4 in a 100x100 lattice. The order parameter 
as a function of 9 from the fit of Fig. 4 (points) and the same quantity using 
the algorithm proposed in jS]. 
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Fig. 6 
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Figure 6: The CP symmetric model presented in the text: numerical (points) 
and exact (line) results. 
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